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Abstract 

The BABAR detector, operating at energies near the T(4S) resonance at the PEP-II asymmetric B 
Factory at SLAC, has collected a sample of 32 million BB pairs by May 2001. A study of time- 
dependent CP-violating asymmetries in events where one neutral B meson is fully reconstructed in 
a final state containing charmonium has resulted in the measurement sin2/3 = 0.59 db 0.14 (stat) ± 
0.05 (syst), which establishes CP violation in the B° meson system at the 4<r level. B lifetime and 
mixing measurements from a sub-sample of 23 million BB pairs are also presented. 
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1 Introduction 



Since its discovery in 1964 in the decays of K® mesons CP violation has been the subject of 
many experiments and the motivation for many theoretical developments in particle physics. The 
three-generation Standard Model accommodates CP violation through the presence of a non-zero 
imaginary phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix 0. The result 
presented here establishes for the first time the existence of CP violation outside the neutral Kaon 
system. 

The primary goal of the BABAR experiment at PEP-II is to perform stringent tests of the 
Standard Model by over-constraining the Unitarity Triangle through CP violation measurements 
(angles a, (3 and 7) and the determination of its sides \V c b\ in semileptonic B decays and 

\Vtd\ in B°B° mixing). 

2 PEP-II 

The PEP-II B Factory Q is an e + e~ colliding beam storage ring complex at SLAC designed to 
produce a luminosity of 3xl0 33 cm _2 s _1 at a center-of-mass energy of 10.58 GeV (T(AS) resonance). 
At the T(AS) resonance B mesons can only be produced as B + B~ or coherent B°B° pairs. The 
time evolution of a coherent B°B° pair is coupled in such a way that the CP or flavor of one B at 
decay time t\ can be described as a function of the other B flavor at its decay time ti and the signed 
time difference At = t\ — ti. The machine has asymmetric energy beams (9.0 GeV electrons on 
3.1 GeV positrons), corresponding to a center-of-mass boost of /?7=0.56. An average separation of 
Pjct ~250 /im between the two B meson decay vertices allows the measurement of time-dependent 
decay rate asymmetries. PEP-II has exceeded its design luminosity by 30% while BABAR, with a 
logging efficiency of >95%, has been accumulating data at daily rates up to 260 pb -1 . 

3 bad ah 

A detailed description of the detector and its performance can be found in Q. The volume within 
the 1.5T BABAR superconducting solenoid contains a five layer silicon strip vertex detector (SVT), 
a central drift chamber with a helium-based gas mixture (DCH), a quartz-bar Cherenkov radiation 
detector (DIRC) and a CsI(Tl) crystal electromagnetic calorimeter (EMC). Two layers of cylin- 
drical resistive plate counters (RPCs) are located between the barrel calorimeter and the magnet 
cryostat. The instrumented flux return (IFR) outside the cryostat is composed of 18 layers of 
radially increasing thickness steel, instrumented with 19 layers of planar RPCs in the barrel and 
18 in the endcaps which provide muon and neutral hadron identification. 

The SVT has a typical single hit resolution of 15 /jm in z and 97% efficiency, while fully re- 
constructed single B decay vertex resolution in z is 50 fim. Charged particle tracking using the 
SVT and DCH achieves a resolution of (ct(pt/pt)) 2 = (0.0015pr) 2 + (0.0045) 2 , where px is the 
transverse momentum in GeV/c. Photons are reconstructed in the EMC, yielding mass resolutions 
of 6.9MeV/c 2 for vr° -> 77 and 10MeV/c 2 for K° s vr° vr°. 

Leptons and hadrons are identified using a combination of measurements from all the BABAR 
components, including energy loss dE/dx in the DCH and in the SVT. Electron identification 
is mainly based on deposited energy and shower characteristics in the EMC, while muons are 
identified in the IFR and confirmed by their minimum ionizing signal in the EMC. Excellent kaon 
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identification in the barrel region is provided by the DIRC, which achieves a K — tt separation of 
>3Aa in the range 0.25-3.0 GeV/c. 



4 Mixing and CP violation with dilepton events 

The mass difference Amj between the two mass eigenstates of the neutral B system (sensitive 
to \Vtd\) is measured by comparing the rate of neutral B meson pairs decaying with the same b 
quark flavor to the rate of decays with the opposite flavor sign. At the T(4S) this time dependent 
asymmetry gives direct access to Am^: 

WAQ - N{B °^ ){At) ~ w^x*') + = cos(Am , . At) (1) 

il ' N(B B )(At) + (N(B B )(At) + N(B B°)(At)) K ' KJ 

where At is the difference between the two neutral B decay times. The simplest way to determine 
the b quark flavor of the decaying neutral B is to use primary leptons as tagging particles. A& md (At) 
can be constructed directly from the rates of "like-sign" (7 + , Z + ) + (l~ , l~) and "unlike-sign" (7 + , l~) 
events. A sample of 100,000 dilepton events has been selected from 23 million BB pairs (charged 
and neutral). Cascade decays background (b — > c — ► I) is suppressed by a neural network using 
kinematic information of the events. The charged B content of the sample is extracted from the 
data together with Am^. In the boost approximation used in this measurement the decay time 
difference is calculated as: At = Az/c < f3j >, where the small flight path of the B mesons 
perpendicular to the z axis is ignored. The Az resolution is extracted from simulation and has 
been validated by comparisons to real data control samples (e.g. J/ip decays). The time dependent 
asymmetry is shown in Fig. |]. The preliminary result is: 

Am d = 0.499 ± 0.010(stat.) ± 0.012(syst.)?ips _1 (2) 

CP and T violation in B°B° mixing can be probed in the same way as in K°K° mixing |1. Its 
magnitude is given by the parameter e: 

Nb(l+,l+)-Nb(l-,l-) ^ 4Re(e) 

MlAtl) = Nb(i+,i+) + m(i-,i-) ~ TTW (3) 

Using our dilepton sample we obtain the preliminary result: 

= ( L2 ± 2.9(stat.) ± 3.6(syst.)) x 10" 3 (4) 

which is the most precise measurement of this quantity to date. The At distribution is also shown 
in Fig. [p. 



5 Time-dependent measurements with exclusively reconstructed 
B events 

High purity B event samples are obtained when the hadronic decay of one B meson is fully recon- 
structed. In such events the kinematics and decay vertex, as well as the flavor or CP content of the 
exclusively reconstructed B are fully determined. The events where in addition the decay vertex 
of the other B meson can be determined are used for lifetime measurements. The B°B° events 
where both the decay vertex and flavor (at decay time) of the other B can be determined are used 
for mixing or CP violation measurements. 
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Figure 1: Left: The time dependent asymmetry between unmixed and mixed dilepton events. 
Right: The asymmetry At- The fit results are superimposed on both asymmetry plots. 

5.1 Exclusive B reconstruction 

B° mesons are fully reconstructed in hadronic modes of (a) definite flavor: B° — > £)W~7r + , p + , 
D^'af, J/^K*°(K*° K+7T-) and (b) known CP content: B° -> J/ipK% tp(2S)K% J/ipK®, 
XdK% J/ipK*°(K*° -» K®ir°). In the following the two samples are referred to as -Bfl av an d Bqp 
respectively. B^ mesons are reconstructed in the hadronic modes B~ — ► D°7r _ , D*°ir~, JfyK~, 
ip(2S)K~ (throughout this paper conjugate modes are implied). 

The selections have been optimized for signal significance, using on-peak, off-peak and simu- 
lated data. Charged particle identification, mass (or mass difference) and vertex constraints are 
used wherever applicable. The signal for each decay mode is identified in the two-dimensional dis- 
tribution of the kinematical variables AE and toes : AE = E* cc — Ef is the difference between the 

B candidate energy and the beam energy and toes = \J 2 — p^ c is the mass of a particle with a 
reconstructed momentum p* ec = ]T\ p* assumed to have the beam energy, as is the case for a true B 
meson. In events with several B candidates only the one with the smallest AE is considered. The 
AE and toes variables have minimal correlation. The resolution in toes is ~3MeV/c 2 , dominated 
by the beam energy spread. The resolution in AE is mode dependent and varies in the range of 
12-40 MeV. For each mode a rectangular signal region is defined by the three standard deviation 
bands in toes (5.27 < toes < 5.29GeV/c 2 ) and AE (mode dependent interval). The composition of 
each sample is determined by fitting the toes distribution for candidates within the signal region in 
AE to the sum of a single Gaussian representing the signal and a background function introduced 
by the ARGUS collaboration @. 

5.2 At calculation and resolution 

Since no stable charged particle emerges from the T(4S) decay point, the production point of the 
B mesons and thus their individual decay times cannot be determined. However the decay time 
difference At between the two is sufficient for the description of a coherent B meson pair (decay 
length difference technique). 

The difference of the proper decay times of the B mesons At =t\ — t% is determined from the 
separation along the boost direction Az =z\ — Z2, including an event-by-event correction for the 
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direction of the B with respect to the z direction in the T{4S) frame. z\ is determined from the 
charged tracks that constitute the exclusively reconstructed B\ candidate. The other B vertex is 
determined by fitting the tracks not belonging to the B\ candidate to a common vertex. Tracks 
from photon conversions are removed. Pairs of tracks compatible with the decay of a long lived K° 
or A are replaced by the parent neutral pseudotrack. To reduce the bias in the forward z direction 
from charm decay products, the track with the largest contribution to the vertex x 2 , if above 6, is 
removed and the fit is iterated until no track fulfills this condition. Knowledge of the beam spot 
location and beam direction is incorporated in the B2 vertex determination through the addition 
of a pseudotrack to its vertex, computed from the B\ vertex and three-momentum, the beam spot 
(with a vertical size of 10 /iin) and the T(4f>) momentum. The Az reconstruction efficiency is 
97%. For 99% of the reconstructed vertices the r.m.s. Az resolution measured in data is 180 yum, 
dominated by the Z2 vertex. An event is accepted if it has converged fits for the two B vertices, an 
error of less than 400 ^m on Az, and a measured \At\ < 20 ps. 

The modelling of the resolution function 1Z is a crucial element of all time-dependent mea- 
surements. In the case of B lifetime measurements, studies both on simulated and real data have 
shown that adding a zero-mean Gaussian distribution and its convolution with a decay exponential 
provides the best compromise between different sources of uncertainties. The width of the Gaussian 
is scaled by the per-event At error, derived from the vertex fits. A zero-mean Gaussian distribu- 
tion with width fixed at 10 ps is used to describe outliers (less than 1% of events with incorrectly 
reconstructed vertices) . 

In the mixing and CP violation measurements the time resolution is described by the sum of 
three Gaussian distributions (core, tail and outliers) with different means. In the CP fit the means 
are modeled to be proportional to the per-event At error, which is correlated with the weight that 
the daughters of long-lived charm particles have in the tag vertex reconstruction. The core and tail 
widths are scaled by the per-event At error. The outlier width is fixed to 8ps. 

5.3 Flavor tagging 

For flavor tagging we exploit information from the incompletely reconstructed other B decay in 
the event. Each event is assigned to one of four hierarchical, mutually exclusive tagging categories 
or excluded from further analysis. The Lepton and Kaon categories contain events with high 
momentum leptons from semileptonic B decays or with kaons whose charge is correlated with the 
flavor of the decaying b quark {e.g. a positive lepton or kaon yields a B° tag). The NT1 and NT2 
categories are based on a neural network algorithm whose tagging power arises primarily from soft 
pions from D* + decays and from recovering unidentified isolated primary leptons. 

The figure of merit for tagging is the effective tagging efficiency Qi = — 2u>j) 2 , where Ej 
is the fraction of events with a reconstructed tag vertex that are assigned to the i th category and 
Wi is the mistag fraction for the same category. The statistical error on sin2/3 is proportional to 
1/s/Q, where Q = YlQi- The efficiencies and mistag fractions for the four tagging categories are 
measured from data and summarized in Table 

5.4 B lifetime and mixing measurements 

The results presented here have been obtained from a sample of approximately 23 million BB pairs 
collected by BABAR between October 1999 and October 2000. Samples of «6000 B° and «6300 B + 
signal events have been selected with background contamination of less than 10%. The results of a 
fit with a Gaussian signal distribution and an ARGUS background function || are superimposed 
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on the m,Es distribution of the final sample in Fig. The B meson lifetimes are extracted from 
unbinned maximum likelihood fits to the At distributions, also shown in Fig. |2| We obtain: 



t b o = 1.546 ± 0.032 (stat) ± 0.022 (syst) ps 
t b + = 1.673 ± 0.032 (stat) ± 0.023 (syst) ps 
t b + /t b o = 1.082 ± 0.026 (stat) ± 0.012 (syst) 




Figure 2: Left: toes distributions of the selected neutral (top) and charged (bottom) B candidates. 
Right: At distribution for the B° (top) and B + (bottom) events within 2a of the B mass in tties 
with superimposed fit results. The single-hatched areas are the background components and the 
cross-hatched areas represent the outlier contributions. 

Neutral B meson pairs, produced as B°B° decay either as B°B° (unmixed) or B° B (B° B ) 
(mixed), allowing to observe B°B° mixing and measure Am^- A sample of ~4500 such events 
where the flavor of the second B at the time of its decay has been tagged has been obtained. The 
At distributions and the mixing oscillation are shown in Fig. ||. From an unbinned maximum 



Tagging Category 




w (%) 


Aw (%) 


Q(%) 


Lepton 


10.9 ±0.3 


8.9 ± 1.3 


0.9 ±2.2 


7.4 ±0.5 


Kaon 


35.8 ±0.5 


17.6 ± 1.0 


-1.9 ± 1.5 


15.0 ±0.9 


NT1 


7.8 ±0.3 


22.0 ±2.1 


5.6 ±3.2 


2.5 ±0.4 


NT2 


13.8 ±0.3 


35.1 ± 1.9 


-5.9 ±2.7 


1.2 ±0.3 


all 


68.4 ±0.7 






26.1 ± 1.2 



Table 1: Average mistag fractions Wi and mistag differences Awi = Wi (B°) - Wi(B°) extracted 
for each tagging category i from the maximum-likelihood fit to the time distribution for the fully- 
reconstructed B° sample (B^+Bcp)- Uncertainties are statistical only. 
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likelihood fit we obtain the preliminary result: 

Am d = 0.519 ± 0.020(stat.) ± 0.016(syst.)7i ps _1 (5) 




Figure 3: Left: At distributions of the selected neutral unmixed (top) and mixed (bottom) events. 
Right: The time dependent asymmetry between unmixed and mixed events. The fit results are 
shown superimposed in all distributions. The dashed curves on the At plots indicate the background 
contributions. 



6 sin2/3 and the observation of CP violation 



The data set of 32 million BB pairs collected between October 1999 and May 2001 has been used 
to fully reconstruct a sample Bqp of neutral B mesons decaying to the J/ipK®, ip(2S)K®, J/ipK®, 
XdKg, and J/ip K*°(K*° — ► K®tt°) final states. The last two modes have been added since our first 
sin2/3 publication 11. There are several other significant changes in the analysis. Improvements 
in track and K® reconstruction efficiency in 2001 data produce a ~30% increase in the yields 
for a given luminosity. Better alignment of the tracking systems in 2001 data and improvements 
in the tagging vertex reconstruction algorithm increase the sensitivity of the measurement by an 
additional 10%. Optimization of the J/ipK® selection increases the purity of this sample. The 
final Bqp sample contains about 640 signal events and, with all the improvements, the statistical 
power of the analysis is almost doubled with respect to that of Ref. ||. A sample Bfi av of 7591 
fully reconstructed B° events in definite flavor eigenstates has been used to determine the tagging 
performance and a sample of 6814 fully reconstructed charged B has been used for validations. 
Details of our samples are shown in Table ^. The samples are shown in Fig. ^. 

We examine each of the events in the Bcp sample for evidence that the other neutral B meson 
decayed as a B° or a B° (flavor tag). The decay-time distributions for events with a B° or a B° 
tag can be expressed in terms of a complex parameter A that depends on both B°B° mixing and on 
the amplitudes describing B° and B decay to a common final state / . The distribution f + (f_ ) 
of the decay rate when the tagging meson is a B°(B°) is given by 



f±(At) = 



-|At|/r s0 



2r B o(l + |A| 2 ) 



IAI 2 1 — IAI 2 
^— - ± 2mA sin (Am B « At) =F ^— - cos (Am B oAt) 



(6) 
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Figure 4: a) Distribution of tties f° r Bcp candidates having a K® in the final state; b) distribution 
of AE for J/ipK® candidates. 



Table 2: Number of tagged events, signal purity and result of fitting for CP asymmetries in the 
full CP sample and in various subsamples, as well as in the -Bfl av and charged B control samples. 
Errors are statistical only. 



Sample 


A^tag 


Purity (%) 


sin2/3 


J/iPK»,iP(2S)KlxdK u s 

J/4>K° ( Vf = +1) 
J/tpK* ,K* -» K° s -k° 


480 
273 
50 


96 
51 
74 


0.56 ± 0.15 
0.70 ± 0.34 
0.82 ± 1.00 


Full CP sample 


803 


80 


0.59 ± 0.14 


J/^Kl 1>(2S)Kl XdK u s only ( Vf 






J/4jK° s {K° s ^7r vr ) 
V>(2£)if° -> vr+vr-) 


316 
64 
67 
33 


98 
94 
98 
97 


0.45 ± 0.18 
0.70 ± 0.50 

0.47 ± 0.42 

o kq _i_ 0.55 


Lepton tags 


74 


100 


0.54 ± 0.29 


Kaon tags 


271 


98 


0.59 ± 0.20 


NT1 tags 


46 


97 


0.67 ± 0.45 


NT2 tags 


89 


95 


0.10 ± 0.74 


5 U tags 


234 


98 


0.50 ± 0.22 


B° tags 


246 


97 


0.61 ± 0.22 


.Bfl av non-CP sample 


7591 


86 


0.02 ± 0.04 


Charged -B non-CP sample 


6814 


86 


0.03 ± 0.04 
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where At = tcp — itag is the time between the two B decays, t b o is the B° lifetime and Amgo 
is the mass difference determined from B°B° mixing || . The first oscillatory term in Eq. || is due 
to interference between direct decay and decay after mixing, and the second term is due to direct 
CP violation. A difference between the B° and B° At distributions or a At asymmetry for either 
flavor tag is evidence for CP violation. 

In the Standard Model A = r\fe~ 2% ^ for charmonium-containing b — > ccs decays, rjf is the CP 
eigenvalue of the state / and /? = arg [ — V cA V* h /V tA V^ ] is an angle of the Unitarity Triangle. Thus, 
the time-dependent CP-violating asymmetry is 

Acp ( At ) = rWj + f~(At| = ~ Vf sin2 ^ sin ( Am ^° At )' ( ? ) 

where r\j = —1 for J/ipK®, i[){2S)K® and x&Ks anc ^ +1 f° r J/^-^-l- Due to the presence of even 
(L=0, 2) and odd (L=l) orbital angular momenta in the J/ifi K*°(K*° — > K®7r°) system, there are 
CP-even and CP-odd contributions to the decay rate. When the angular information in the decay 
is ignored, the measured CP asymmetry in J/ip K*° is reduced by a dilution factor D± = 1 — 2Pj_, 
where R± is the fraction of the L=l component. We have measured R± = (16 ± 3.5)% || which, 
after acceptance corrections, leads to an effective rjt = 0.65 ± 0.07 for the J/tpK*° mode. 

The Bqp and Bfi av samples are used together in the unbinned maximum likelihood fit for the 
extraction of sin2/3. A total of 45 parameters are varied in the fit, including sin2/3 (1), the average 
mistag fraction w and the difference Aw between B° and B° mistags for each tagging category 

(8) , parameters for the signal At resolution (16), and parameters for background time dependence 

(9) , At resolution (3) and mistag fractions (8). The determination of the mistag fractions and 
signal At resolution function is dominated by the large -Bfl av sample. Background parameters are 
governed by events with w-es < 5.27GeV/c 2 . As a result, the largest correlation between sin2/3 
and any linear combination of the other free parameters is only 0.13. We fix t b o = 1.548 ps and 
Am B o = 0.472 ftps^ 1 ||. The value of sin2/? and the CP asymmetry in the At distribution 
were hidden (blind analysis) following our previous publication ||, until the event selection was 
optimized and all other aspects of the present analysis were complete. 

Figure || shows the At distributions and Aqp as a function of At overlaid with the likelihood 
fit result for the rjf = —1 and r\$ = +1 samples. The probability of obtaining a lower likelihood, 
evaluated with a parameterized simulation of a large number of data-sized experiments, is 27%. 
Our result is: 

sin2/? = 0.59 ± 0.14 (stat) ± 0.05 (syst). 

Repeating the fit with all parameters except sin2/3 fixed to their values at the global maximum 
likelihood, we attribute a total contribution in quadrature of 0.02 to the error on sin2/3 due to the 
combined statistical uncertainties in mistag fractions, At resolution and background parameters. 
The dominant sources of systematic error are the parameterization of the At resolution function 
(0.03), due in part to residual uncertainties in SVT alignment, possible differences in the mistag 
fractions between the Bcp and -Bfl av samples (0.03), and uncertainties in the level, composition, 
and CP asymmetry of the background in the selected CP events (0.02). The systematic errors 
from uncertainties in Am B o and t b o and from the parameterization of the background in the -Bfl av 
sample are small; an increase of 0.02 ?ips -1 in the value for Am B o decreases sin2/3 by 0.015. 

The large sample of reconstructed events allows a number of consistency checks, including 
separation of the data by decay mode, tagging category and i?tag flavor. The results of fits to these 
subsamples are shown in Table [2|. Table |2| also shows results of fits to the samples of non-CP decay 
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Figure 5: Number of rjf = —1 candidates (J/ipK®, Tp(2S)K®, and Xc\K®) in the signal region a) 
with a B° tag N B o and b) with a B° tag Ng , and c) the asymmetry (N B o — N^ )/(N B o +N^ ), as 
functions of At. The solid curves represent the result of the combined fit to all selected CP events; 
the shaded regions represent the background contributions. Figures d)-f ) contain the corresponding 
information for the rjf = +1 mode (J/ipK®). The likelihood is normalized to the total number of 
B° and B° tags. The value of sin2/3 is independent of the individual normalizations and therefore 
of the difference between the number of B° and B° tags. 



modes, where no statistically significant asymmetry is found. Performing the current analysis 
on the previously published data sample and decay modes yields a value of sin2/3 =0.32±0.18, 
consistent with the published value |J. For only these decay modes, the year 2001 data yield sin2/3 
=0.83±0.23, consistent with the 1999-2000 results at the 1.8cr level; for the J/ip K° S {K° S vr+vr-) 
channel the consistency is at the lAa level. 

If |A| is allowed to float in the fit to the rjf = — 1 sample, which has high purity and requires 
minimal assumptions on the effect of backgrounds, the value obtained is |A| = 0.93 ± 0.09 (stat) ± 
0.03 (syst). The sources of the systematic error in this measurement are the same as in the sin2/3 
analysis. In this fit, the coefficient of the sin (Am^oAf) term in Eq. ^ is measured to be 0.56 ± 0.15 
(stat) in agreement with Table [2|. 

The measurement of sin2/3 = 0.59 ± 0.14 (stat) ± 0.05 (syst) establishes CP violation in the 
B° meson system at the 4.1cj level. This significance is computed from the sum in quadrature of 
the statistical and additive systematic errors. The probability of obtaining this value or higher 
in the absence of CP violation is less than 3 x 10~ 5 . This direct measurement is consistent with 
the range implied by measurements and theoretical estimates of the magnitudes of CKM matrix 
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7 Conclusions and prospects 

37 years after the discovery of CP violation in the Kaon system, BABAR has established CP violation 
in the B system with the measurement: 

sin2/? = 0.59 ± 0.14 (stat) ± 0.05 (syst). 

By the summer of 2002 BABAR will have a data sample of more than 100 million BB pairs, bringing 
the precision on sin2/3 to less than 0.1 and allowing searches for other manifestations of CP violation 
in the B system. 
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